Neutron magnetic form factor in strongly correlated materials 
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We introduce a formalism to compute the neutron magnetic form factor FA/(q) within a first- 
principles Density Functional Theory (DFT) + Dynamical Mean Field Theory (DMFT). The ap- 
proach treats spin and orbital interactions on the same footing and reduces to earlier methods in 
the fully localized or the fully itinerant limit. We test the method on various actinides of current 
interest NpCoGas , PuSb and PuCoGas , and we show that PuCoGas is in mixed valent state, which 
naturally explains the measured magnetic form factor. 

PACS numbers: 71.27.+a, 74.20.Mn, 75.25.-j 



Compounds including elements from the actinide series 
provide a beautiful illustration of the challenges posed by 
correlated materials. The 5/ electrons in these systems 
display simultaneously itinerant (i.e. band-like) and lo- 
calized (atomic-like) properties. Describing the impact 
of this wave-particle duality on different physical observ- 
ables, measured using different spectroscopic probes, is 
an outstanding theoretical challenge. 

Neutron scatteringllj is a time-honored probe to inves- 
tigate the dynamics of the magnetic degrees of freedom. 
It probes the dynamic susceptibility, describing the spa- 
tial and temporal distribution of magnetic fluctuations. 
In the itinerant limit, it can be modeled in terms of a 
particle hole continuum of quasiparticles, while in the lo- 
calized limit it can be describe in terms of propagating 
spin waves. It is generally accepted that in many mate- 
rials neither a fully itinerant nor a fully localized picture 
is adequate and some combination of both is required to 
model the dynamics of the spin fluctuations as in the du- 
ality model of Ref. 12- 

The intensity in the magnetic Bragg peaks can be used 
to obtain a real picture of the magnetization inside the 
unit cell. This can be done even for materials that do 
not exhibit magnetic long range order, by applying an 
external magnetic fleld. Classical techniques can handle 
a fully itinerant or a fully localized picture [3|. However 
these approaches are not sufficient for many compounds 
of considerable scientific interest. It has been known for 
a while that intermediate valence rare-earth semiconduc- 
tors show puzzling magnetic properties that can be ex- 
plained only by a theory which explicitly considers the 
spatial extend of the magnetic excitations [J| . Similarly 
only magnetic orbitals of strong covalent nature can cor- 
rectly account for the neutron intensity in the cuprates 
[5|. A theory able to describe the magnetic form fac- 
tor for partly itinerant systems from first principles is 
needed. 

Important recent experiments of Hiess et. al. de- 
termined the magnetic field induced form factor of 
PuCoGas, a material which superconducts at the remark- 



ably high transition temperature Tc ~ 18.5 K, a record 
in the heavy- fermion family [6| . The degree of itinerancy 
of the f electrons is the subject of active debate and has 
important consequences for the mechanism of supercon- 
ductivity. Neither the localized nor the itinerant model 
of the neutron form factors fits the data well, providing 
strong motivation for our theoretical developments. 

In this letter we develop a method to compute 
the form factor for magnetic neutron scattering within 
DFT-HDMFT 0. We test the method on several ac- 
tinide materials. The FuCoGas induced magnetic form 
factor is consistent with correlated mixed valent nature of 
the material, where both the 5/^ and 5/^ configuration 
are important. This is reminiscent of the mixed valent 
nature of elemental plutonium [8| . 

The magnetic form factor FM(q) is defined by 



rM(q) = - 



1 



2/ii3 



MT(q) 



(1) 



where MT(q) = MjP'"(q)+M^''^(q) is the Fourier trans- 
form of the transverse component of the magnetization 
density q x (M(r) x qj, fiB is the Bohr magneton and 
q is the scattering wave vector at the Bragg peak. To 
avoid ambiguity in definition of magnetization [9[, we 
express the form factor in terms of the Fourier trans- 
form of the current density J(q) — J dr e~''''''J(r). The 
current and the transverse magnetization are related by 
MT(q) = ^q X J(q)/(7^. The current has two contribu- 
tions, the spin part Jspin(r) and the orbital part Jorb(r)- 
Expressing the definition of Jorb(r) and Jspin(r) in terms 
of fleld operators ^^(r), we flnd for the form factor the 
following expression 



FM(q)-^E/dr'^"">-^ 



*l(i-)qx 



ffss' X q + ^^s/V 



(2) 



*s'(r) 



where c? is the vector of the Pauli matrices and s, s' are 
the spin indexes. For a more detailed derivation see the 



on-line supplementary material [10| . It is useful to notice 
that the limit limg_^oFM(q) can be well defined, but it 
is subtle [IjJ. However the form factor FM(q) is mea- 
sured only at finite q values and hence it is free from 
ambiguities. 

rM(q) is measured in polarized-neutron diffraction ex- 
periments directly through the flipping ratio technique. 
In this method an external magnetic field B is applied 
to the sample, and the ratio R = {da/dfl)_^_ / {d<j/<Kl)_ 
between the cross section for neutrons polarized par- 
allel and anti-parallel to B is measured. In a cen- 
trosymmetric crystal structure with collinear magnetic 
moments and q _L B the flipping ratio R satisfies 
I Sr'^'^ ) = 7?'o-FA/(q)/& where b is the known nuclear 

scattering amplitude, Fm{(\) the component of the mag- 
netic structure factor parallel to B, 7 = 1.9132 and 
To — fie^ jmc? is the classical electron radius. More gen- 
eral formulas which relate the form factor to the flip- 
ping ratio for other crystal structures and experimen- 
tal setups are given in Ref. [l2- For localized electrons 
the form factor is commonly fitted to the following ra- 
dial dependence Fm(<7) = -2^((jo(9)) + C'2(J2('7))), 
where {jkil)) stands for the spatial average over the 
atomic wave function of the magnetic atom (which is 
usually solved in the isolation). This should be under- 
stood in the so called dipole approximation. The expo- 
nent e"*'''' is expanded around the center of the atom as 
e"'''' « jo{qr) - i{q- r){jo{qr) + J2{qr)), where jk{qr) 
are spherical Bessel functions of order k. Within this ap- 
proximation, the form factor is greatly simplified and in 
the common experimental set up (q _L B, B = Bz,), it 
reduces to 



FAliq) 



Szjoiqr) + -Iz {jo{qr) +J2{qr)} ) . (3) 



Here r is the distance from the magnetic atom, and 
(• • • ) stands for the spatial and temporal average. The 
first and the second term in Eq. (j3]) come from the 
spin and the orbital contribution, respectively. The 
comparison of the above expansion with Eq. ([3]) shows 
that ^ = — /is(2s2 + Iz) and /iC2 = —fJ-silz), hence 
C2 = tJ-L/if-L + fJ-s)- Clearly the ratio C2, which is 
given by the shape of the form factor, uniquely deter- 
mines the size of the orbital and spin component within 
the dipole approximation. Even so, caution is necessary 
in interpreting experiments with Eq. ([3]), since a priori 
the magnitude of higher order terms beyond the dipole 
approximation is not known 



13|,llJ 



To compute the form factor within DFT+DMFT, we 
apply a small magnetic field B = Bz to induce a finite 
magnetic moment. We solve the DMFT problem in the 
presence of magnetic field, and evaluate the correlation 
function Eq. 1^. When expressed in the Kohn-Sham 



basis, Eq. ([2]) takes the form 

T? fr^\ ^ \^ „DMFT 



X / dre^"i>ki(i':s) qx 

/unit cell 



1 



-ass' X q+5syV 



V'kj(r,s')7 
(4) 



where ipki(r,s) are the Kohn-Sham orbitals, i runs over 
the Kohn-Sham bands, and k over the first Brillouin 



zone. The "DMFT density matrix" 



^DMFT 



is expressed 



in terms of the DMFT Green function Gij (k, w) in the 



solid 



,DMFT 



J_ 

2iri 



Jdoj {G*^{Koj) 



G,,(k,L.))/(w) 



where /(w) is the Fermi function. The form factor is 
thus expressed in terms of the one particle correlation 
function, which is easily accessible within DMFT. More- 
over, the spatial integral is local and runs over one unit 
cell, which makes local DMFT approximation particu- 
larly suitable for this problem. We implemented Eq. (|4]) 
within the recent realization of DFT-f-DMFT [l5| based 
on Linear Augmented Plane Wave (LAPW) basis set 
as implemented in the full potential electronic structure 
code Wien2k [16[. The explicit formulas for the form 
factor evaluation within this basis set, as well as detail 
derivation of Eq. (|4]) are given in the on line material 
[lOj. To solve the impurity problem in the presence of 
magnetic field, we used the Non-Crossing Approxima- 
tion [15| . Our calculations show small anisotropic correc- 
tions to dipole approximation for the materials studied 
here, suggesting that the dipole approximation is a good 
approximation for these compounds. For comparison, we 
also compute the form factor within Local Spin Density 
Approximation (LSDA) as first discussed in Ref. |l7|. In 
practice we evaluate the mean value of Eg. (EI inside the 
atomic sphere following the lines of Ref. |l8|. We perform 
the LSDA calculation in the presence of external mag- 
netic field, as implemented in Wien2K [19|. 

In Fig.mA) we compare theoretical DFT-hDMFT and 
LSDA form factors with experiments on NpCoGas in the 
paramagnetic state [y] . Our DFT-I-DMFT form factor is 
in excellent agreement with experiment, while the LSDA 
dramatically fails in this material. The LSDA form factor 
shows a minimum at finite wave vector q. Such a large 
minimum can be explained by C2 ~ —9.5; this occurs 
since fiL and fig almost cancel, but l/i^l < \fis\- An un- 
derestimation of the orbital moment is typical of LSDA. 
Within DFT-I-DMFT the atomic degrees of freedom are 
treated exactly by the exact diagonalization of the atomic 
5/-shcll in the presence of magnetic field. This ensures 
that Hund's rule coupling is properly treated, leading to 
anti-parallel fiL and fis, but |/iL| > |//s|, hence C2 > 0. 
For NpCoGa5 we determine the value of the coefficient 
C2 = 2.16. This value is consistent with localized 5/- 
electrons in the configuration 5/*, in agreement with 
Mossbauer spectroscopy [20[ and neutron diffraction ex- 
periments [6|,|21|. At the same time NMR 22[ and inelas- 
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FIG. 1: (Color on line) A) panel: magnetic form factor 
for NpCoGas. Red dots are experimental data reproduced 
from II. The blue curve with squares is the DFT+DMFT 
calculation, and the black curve with triangles is the LSDA 
calculation. The DFT+DMFT form factor agrees with ex- 
periment with value of the Pearson correlation coefficient 
i?PMCC = 0.95. B) panel: spectral function Ajmj {^) for Np 
f-electrons. Blue curves correspond to the j = 5/2 multiplet 
and green curves to the j = 7/2 multiplet. The experimental 
and theoretical DMFT temperature is T = 52 K. 
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pressed with the consequent absence of a quasiparticle 
multiplet structure in the spectral function [27[ . This re- 
sult is consistent with neutron diffraction data: the form 
factor curve has a characteristic maximum at finite q, 
feature typical of a pure /^ configuration state for the 
Pu atom [2JI . The LSDA calculation underestimates the 
orbital moment and finds a negative C2 coefficient. Our 
DFT+DMFT calculation reproduces the /-electrons oc- 
cupation value (uf) ~ 5.0 of the previous experimental 
and theoretical works 24|, |27[ and indeed it is in good 
agreement with the measured data, see (see Fig. dJ^.). 
In particular we find that there is a large cancellation 
between orbital and spin moment with fJ-s/ fJ-L = —0.74 
and C2 = 3.92. We now turn to PuCoGa.g. Photo 
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FIG. 3: (Color on line) A) panel : magnetic form factor for 
FuCoGas . The blue curve with squares corresponds to the full 
DFT+DMFT calculation, the black curve with triangles to 
the LDA calculation. Red dots are experimental data [a| The 
DFT+DMFT curve agrees with experiment with a Pearson 
correlation coefficient i?pMCC = 0.70 . B) panel: spectral 
function Ajmj (ijj) for Pu f-electrons. The color legend is the 
same as in Fig. [l] 



FIG. 2: (Color on line) A) panel: magnetic form factor for 
PuSb at T = 20 K. Red dots are experimental data [2J], the 
blue curve with square is the DFT+DMFT calculation and 
the black curve with triangles the LSDA calculation. The 
DFT+DMFT curve agrees with experiment with a Pearson 
Correlation coefficient i?pMCC = 0.95. B) panel: Spectral 
function for PuSb. The color legend is the same as in Fig. [l] 



tic neutron scattering [23[ suggest that NpCoGas shows 
also itinerant aspects of the 5 /-electrons. A signature 
of this moderate delocalization is also apparent in our 
calculated spectral function a.t T = 52 K displayed in 
Fig. lUB). A small quasiparticle peak is formed at the 
Fermi level, a signature of electron itinerancy at low en- 
ergy. Next we compute the form factor for PuSb in the 
ferromagnetic state. PuSb is a metal [25|, which orders 
antiferromagnetically below Tjv = 85 K and becomes 
a ferromagnct at T = 67 i^ [26|. Theoretically it has 
been showed that in PuSb valence fluctuations are sup- 



emission spectra show the formation of a quasiparticle 
peak at the Fermi level, however there is a large discrep- 
anc y in the peak height between different measurements 
28l |29| . First magnetic susceptibility measurements sug- 



gested that 5/-electrons behave as unquenched local mo- 
ments until they enter in the superconducting state |30| . 
In turn neutron scattering shows a temperature inde- 
pendent magnetic susceptibility, implying the absence of 
magnetic moments such as in ^-Pu [6|, |31|. Electronic 
structure calculations qualitatively support the picture 
of delocalized 5/-states, however they predict a Pu ion 
close to magnetic order and a form factor shape not ob- 
served in experiments [32|, |33[ . Since our understanding 
of superconductivity in PuCoGas depends on the itiner- 
ant or localized nature of correlated electrons [3^, fur- 
ther theoretical and experimental investigations are com- 
peUing. Within our DFT+DMFT calculation we find 
that a quasiparticle peak appears at the Fermi level, see 
Fig. [DJB). These results are consistent with a specific 



heat coefficient 7 ~ 70 mJ/{K'^mol), which compares 
weU with experiments (30[, and go beyond the pioneer 
DFT+DMFT calculations, solved within the T-matrix 
and fluctuating exchange technique [35|. Together with 
a quasiparticle peak, a mixed valent state forms, where 
the 5/-electrons have a finite probability to be both in 
the configuration state /^ and /^. Our theoretical pre- 
diction for the 5/^ occupation probability is Pje = 0.26, 
corresponding to (?i/) ~ 5.26 and a coefficient C2 = 2.35. 
We plot the corresponding form factor curve in Fig. ^A) 
together with the form factor obtained from the LSDA 
calculation. As for the previous materials, LSDA under- 
estimates the orbital moment and it obtains a negative 
C2 coefficient that is inconsistent with experimental data. 
The DFT-hDMFT form factor with C2 = 2.35 well de- 
scribes the neutrons data and it accounts also for the 
magnetic susceptibility ( see the supplementary mate- 
rial [10|). The value of C2 = 2.35 is naturally explained 
by the mixed valence picture obtained theoretically for 
PuCoGas. For a free Pu'^"'" ion solved in the intermedi- 
ate coupling C2 = 3.83, hence ^-hh^-s = —1.83 jsl; As a 
mixture of the configuration /^ is included in the many 
body ground state, the ratio ^L/fJ-s becomes more nega- 
tive and therefore C2 decreases. As pointed out in Ref. la 
the FM{q) shape is very different from the one expected 
for a pure 5/^ configuration of an isolated Pu ion, as for 
example is found in PuSb, see Fig.[2t^a). At the same time 
it is very different from the LSDA prediction. Hence, the 
magnetic properties of PuCoGa.s are not captured either 
by a free moment picture or by an itinerant picture. 

In conclusion in this letter we presented a new ap- 
proach to compute the neutron magnetic form factor. 
The LSDA treatment fails to reproduce the correct form 
factor since the exchange energy is orbital-independent 
and therefore Hund's rules are not respected. On the con- 
trary DFT-I-DMFT includes the atomic physics needed 
to describe strongly correlated systems. Application of 
DFT-KDMFT to PuCoGas suggest an explanation of the 
results of Ref. |6| in terms of a mixed valence picture 
where the ground state of Pu fluctuates between two 
distinct configurations: p and /^. We indeed checked 
that this picture accounts for the values of the specific 
heat and susceptibility as well as for the shape of photo- 
emission spectra. We find a close similarity between the 
DFT4-DMFT valence histogram of PuCoGas and (5-Pu, 
suggesting a close analogy of the local physics in these 
two materials; the magnetic form factor of PuCoGas 
would then be very similar to that of J— Pu, for which 
experiments are notoriously difficult. Finally, mixed va- 
lence is an attractive mechanism for pairing in heavy 
fermions [36[, which could account for the high temper- 
ature superconductivity in PuCoGas. Acknowledgment: 
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Supplementary Material 



FORM FACTOR IN SECOND QUANTIZATION 



The neutron magnetic form factor can be evaluated in terms of current J^f) by the following formula 

imp 1 



FMiq) 



-gx (J(g)), 



eh q^ 
where me, e are respectively the electron mass and the electron charge (e < 0), if is the scattering wave vector and 



(SI) 



Jiq) = J dre-*''-'^(Jorb(r) + Jspin(r)) . 

The orbital and spin currents in terms of the field operators ^s(r') are given by 

eh 



Joi-h{r) 



2im 
ih 



■^ [vI't(f)(VM/,(f)) - (Vvl/t(r))*,(f) 



4in(r) = ^^[*t(r)(V*,,(r)) + (V*l(r))*.,(r) 



(S2) 

(S3) 
(S4) 



Note that the Fourier transform of the spin current Js-pin{q) greatly simplifies: in fact integration by parts leads to 
cancellation of the derivatives of the field operator and we obtain Js-pin{<l) = —ehi/{2me) X^^s' '^ss' Pss'io) x Q- 
Inserting the expression for the currents into Eq. (jSip . we get 



^^M ^^T.J dre-^'^'i'lir) Q 



-ass' X q + Sss'"^ 



*.'(^) 



(S5) 



Now we express the field operator in terms of a complete set of one particle wave functions, such as Kohn-Sham 
orbitals 



*s(r)= J2 ^k^{r,s)cj:^ 



iGbnd,kGlBZ 



to get 



^'(^^^E"£rl,/ 



k,ij,ss' 



dfe "^'Vj-C^, s) g X 



-ass' X q + Sss'^ 



i^kAr,s') 



(S6) 



(S7) 



Here n?^^^'^ is the DMFT density matrix expressed in the Kohn-Sham base, k runs over the first Brillouin zone only, 
and i,j run over Kohn-Sham bands. Because qis reciprocal vector, the integration over space is performed only over 
one unit cell (denoted by cell). When the spin-orbit coupling is large, bands do not have the spin index, because there 
is mixing between both spins species, hence we need to double the number of bands. 
In the LAPW basis, there are two contributions to the above equation: within muffin-tin and in the interstitial region. 



Within Muffln-Tin 



Inside the Muffin- Tin the KS orbitals are expressed in terms of coefficients A - and LAPW basis functions 

i,Ks 



Xfe+/?(^'«) 



as 



V'felr, s) = ^ a1^^xuk(^^ *) = E 4ks 41 <*('^*) yL{n)xs , 



(S8) 



K 



Ks 



where K are reciprocal lattice vectors, t marks the atom type, u°*(r), uj*{r), uf*{r) are the radial solutions of the 
Dirac equation, its energy derivative, and optional local orbitals; YL{ft) are spherical harmonics with L labeling the 
angular quantum numbers l,m . 
For shorter notation, we define a new type of density matrix in the muffin-tin subspace 



^DMFT 
k.t.LsH 



sk.L's'k' Z-^ I ^—^ t.Ks LK I k,ij I ^—^ J,Ii s' L'K' 

-i-j \ K I \k' 

and we express the form factor inside the mufhn-tin in terms of this density matrix 

"1 



F, 



M 



w = ;Jr E 



q - 

k,t,LsK,,L' s' K' 



k,t,LsK,L's'K,' ^ ' -1^ A* I 



qx {uss' X q) + Sss'qx V 



(S9) 



K Yl' Xs')t (SfO) 



Here integration runs over muffin-tin sphere t. 

From the above expression it is not obvious that the (7 — > is well behaved. However, we can add any constant to 
exponent e*''" in the orbital part of the expression, because it vanishes due to symmetry, ft is therefore possible to 
write an alternative expression 



pA/T/-A _ 1 V^ „DMFT /„,K V 

which clearly is well behaved in the q ~ limit 



L Xs 



ie-'9~<f X (a,,, xq) + <5,,,(e-''9"- l)qx V 
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Interstitial Region 



In the interstitial region the KS solution is 



T^ki\ ' J / J i.Ks , /T/ ,, K 



{k+K)r^ 
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We again define a corresponding density matrix 
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and express the form factor by 
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^K'-K-q~Yl 
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(SI4) 



J Vceii \K' -K -q\ 
The last term comes from the difference of the integral over the entire unit cell and inside all muffin-tin spheres. 



DIPOLE APPROXIMATION 



In the dipole approximation, we approximate the Fourier exponent with 



and obtain 






q - 

k,ij,ss' 



'■ Jo{qr) - J((f • r){jo{qr) + J2{qr)) 



dr[jo{qr) - i{q- r){jo{qr) + J2{qr))] V'L(^s) 9 x 



ass' X q + Sss'V 



(S15) 



^f:^{r,s'){SlQ) 



Due to parity selection rules, only the following two terms arc nonzero 



k.ij.ss' 



X / dri/jUr^s) 



cell 



ki 



Joiqr)-qx {ass' x g) - iSss' [jo{qr) +J2{qr)]iq- f)(g x V) 



V'fc,(^,s')- 



Inside the expression Eq. (|817p for form factor, we can use 



i{q-r){qx V) = 2~^di^^^'^^^^ ^"^ ^ l^^iq ' ^{q ^ V) + -i(q • V)(g x f) . 



The following term vanishes 



k,ij 



Hence we can use 



1 ^ 1 ^ i 

-i{q- r}{qx V) + -i{q- V)(g x r) = -qx [qx (f x V)] 



(S17) 



(S18) 



(S19) 



and f X V = T-/ to write 

h 

i{q- r){qx V) = —qx (fx cf) ; 
Finally we find the expression for the form factor in the dipole approximation 



q - 

k^ij^ss' 



^DMFT 
k,ij 



X / dfi'Uf,s) 



-qx {gss' X q}ja{qr)+Sss'^[iqx I) x g] [jo M + J2 M] 



V%-(^-s')- 



(S20) 



(S21) 



Inside the muffin-tin the Kohn Sham orbitals are expressed in terms of a basis in which I and s are diagonal. Therefore 
the form factor in the dipole approximation, inside the muffin-tin, reads 



FM{q) = E 



^^DMFT 

k,t,LsK,LsK' 



qx{ezXq) 



k^t.sL^K.hi' 

where we took the magnetic field in z direction. 



SzK\]o{qr)\uf ) -f -h (ufljoiqr) + J2iqr)\uf ) 



(S22) 



NUMERICAL EVALUATION 



For numerical evaluation we split the form factor expression into the dipole part and the correction to the dipole 
approximation AFm- The correction is 

k,ij,ss' 

Inside the muffin-tin sphere, this expression takes the form 



^FM{q-) = ^ _ E -fZL... {-^x (-^^ X 'r)rls{uiY,\ [e-'''-Mqr)] \u,Y„) 



(S23) 



k.t.s.LK.L' K 



{uiYlI [e 



l + i 



i{q- ^ {jo{qr) + J2{qr))] {qx V)!^,.^^')} 



Wc compute in advance the following quantities 

aLL'ir) = {YL\e-'^" - 1\Yl') (S24) 

bLL'ir) = {YL\-iq-f]YL') (S25) 

SLL'ir) = {YL\qxV\YL') (S26) 

dLL'ir) = {YL\qxep\YL') (S27) 
and obtain the following one dimensional integral over the radial distance in the mufhn-tin sphere 

^PMiq) = ^ E ^IZI.L's.'t {9 X (^^- X q)\al j drr^uf{r)uf: {r)[aLL.{r) + 6ll- (1 - JoM)] (S28) 

k,t,s,LK.,L' K.' 



+ /'drr2<(r)^^^E{«^^"('') " ^""W Mv) + 32{qr)]} d^^' u{r) 



dr 

L" 

+ I drr^u^{r)uf'{r)"^{aLL"{r)-bLL"{r)[jo{qr)+j2{qr)]}cL"L'{r)\ . 
To check numerical accuracy, one could check the accuracy of the following sum 

J2bLL"CL"L' = {YL\-l{q-f)iqX V)\Yl') = 
L" 

-i{YL\{qxr){q-V)-tqx {q X r)\YL') = 
-i{YL\iqx r){q- V)\Yl.) - qx {qx bDUlL' ; 
since for L — L' the first term is zero, we have 

J2bLL"CL"L = -qx{qX e,)h (S29) 

L" 

To derive the above equation, it is useful to know the following property of the spherical harmonics 

^y,,™(0,(^) = ^^{l-m){l + m+l)e-'^Yi^^+^{e,cl,) ~ i ^G + H(^ - "^ + l)e**>l,m-i(e, 0) (S30) 

which follows from 

^l_a;2i.p (3,) = _ Ip (2;) + i(; + ^)(/ _ ^ + l)P,,„_i(a;) . (S31) 

ax 2 2 

MAGNETIC SUSCEPTIBILITY FOR PuCoGas 

From the calculation of the magnetic form factor Fmiq) we can extract the magnetic susceptibility as x = f^/B, 
where fi = —2^bFm{Q) is the magnetic moment of the Pu atom and B is the applied magnetic field. The results are 
summarized in table HI 

TABLE L Magnetic susceptibility for PuCoGas obtained with the NCA impurity solver. 



T (Kelvin) 



12.5 
25 
50 



X (/ia/Tesla) 



12 X 10"'* 

13 X 10"" 
15 X 10"" 



The computed magnetic susceptibility compares well with the values obtained by the neutron experiment [ A. Hiess et. 
ai, Phys. Rev. Lett. 100, 076403 (2008)]. DFT+DMFT obtains a larger susceptibility that the measured one, which 
is comprehensible since it is well known that the Non Crossing Approximation (NCA) impurity solver underestimates 
the Kondo temperature. 

To solve the impurity problem in presence of an applied magnetic field we gave special attention to the off-diagonal 
terms in the impurity hybridization strength A, which give a significant contribution to the form factor. 



